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Porous medium (PM) engine was a new type engine based on the technique of combustion in porous 
medium, which can realize homogeneous and stable combustion. In this paper, the combustion and 
working processes of a specific PM engine were simulated by a two-zone model considering the influ¬ 
ences of the mass distribution, heat transfer from the cylinder wall, mass exchange between zones and 
the heat transfer in porous medium. Influences of operating parameters, e.g. intake temperature and 
pressure, compression ratio, the excess air ratio on the performance of the PM engine were discussed. 
It was found out that the porous medium, acting as a heat recuperator, can significantly enhance the 
evaporation of liquid fuel and preheat the mixture, which promotes the ignition and combustion in 
the cylinder; and that the initial PM temperature and the compression ratio are critical factors controlling 
the compression ignition of the mixture. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

Homogeneous charge and compression ignition (HCCI) is one of 
the most attractive combustion modes for automotive engines to 
achieve both high efficiency and low NO x and particulate emis¬ 
sions. Despite many advantageous features of the HCCI combus¬ 
tion, it still fronts some challenges including higher HC and CO 
emissions and expanding operating range, As no external means 
are used to ignite the mixture, it is especially difficult to control 
the combustion processes inside the engine [1]. In this respect, 
the new concept of the so called porous medium (PM) Engine, 
based on the regenerative or super-adiabatic combustion in a por¬ 
ous medium, offers new ways in dealing with these problems. Re¬ 
cently, the PM engine has received attention from numerous 
researchers because of its potential for producing homogeneous 
mixtures and reducing NO* and soot emissions [2,3]. 

With HCCI, the start of combustion is mainly controlled by the 
chemical kinetics, which depends on the charge composition and 
the pressure and temperature histories of the reactants during 
the compression process. Thus a detailed chemistry mechanism 
is essential for modeling HCCI combustion [4,5]. Various types of 
models have been developed to simulate HCCI combustion includ¬ 
ing single-zone, multi-zone and multi-dimensional models. 

Single-zone models [6] treat the entire cylinder as a homoge¬ 
neous reactor with uniform temperature, pressure and composi¬ 
tions. Because single-zone models represent the highest 
temperature regions in the cylinder, the burn point can be well 
predicted within different operating ranges. However, the inhomo¬ 
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geneity of temperature and compositions within the cylinder are 
not taken into account, so peak pressure, burn duration and the 
amount of emissions could not be accurately predicted. Multi-zone 
models [7,8] divide the cylinder into a number of zones with differ¬ 
ent temperatures and compositions. With the interaction between 
zones and detailed chemical kinetics mechanism being considered, 
the combustion process is reproduced. Multi-zone models can, 
with moderate computation costs, provide results in better agree¬ 
ment with experiments than single-zone models due to consider 
temperature and compositions inhomogeneity in the cylinder. Be¬ 
cause of these advantages the multi-zone models have found wide 
applications in the engineering. Multi-dimensional models [9,10] 
combine detailed chemical kinetics and engine computational fluid 
dynamics (CFD) to simulate HCCI combustion. However, a coupled 
CFD and detailed chemistry simulation requires substantial mem¬ 
ory and CPU time which might be impractical for current computer 
capabilities. 

Development of porous burners has been encouraged by lower 
emission standards in recent years. As a burner of liquid fuel, por¬ 
ous medium is also an efficient evaporator. Combination of large 
heat capacity of the porous material, large specific surface area 
with excellent heat transfer in PM volume makes the fuel droplets 
vaporization very fast and complete. The collision of the fuel drop¬ 
lets within the complex porous structure contributes to very effec¬ 
tive mixture formation and homogenization, which ensures a 
homogeneous combustion in the PM volume. Liquid vaporizing 
combustion in porous medium burner has fine flame stability 
and characteristic of low emission, which has been proved by 
experimental and numerical research [11-13]. Sumrerng and 
Narongsak [11] built a simple porous burner system for burning 
kerosene. The evaporation mechanism and the combustion 
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characteristics occurring in the burner were studied by measuring 
thermal structures in terms of temperature profiles along the cen¬ 
terline of the porous burner. Echigo and Martynenco [12] studied 
the phenomenon of self-sustaining combustion of a gaseous mix¬ 
ture in inert highly porous medium with prior vaporization of li¬ 
quid droplets numerically. Combustion and several heat transfer 
modes as well as droplet behavior in porous medium were inves¬ 
tigated combinedly. Kayal and Chakravarthy [13] developed a 
one-dimensional heat transfer model to analyze the combustion 
of liquid fuel (kerosene) inside an inert porous medium under stea¬ 
dy state conditions using a single step global reaction mechanism. 

Based on these previous researches, a new concept of controlla¬ 
ble combustion in porous media for internal combustion engine 
has been suggested and developed [14-16]. Hanamura et al., [14] 
designed a reciprocating heat engine, similar to a Stirling engines, 
based on the technique of super-adiabatic combustion in porous 
media. One-dimensional numerical simulations showed that the 
thermal efficiency of the engine reached to 57.5% with a compres¬ 
sion ratio between two and three, which is much lower than those 
of conventional Otto and diesel cycles. Weclas [15] suggested a 
strategy for development of intelligent combustion systems for IC 
engines, whose essence involved a new concept for mixture forma¬ 
tion and homogeneous combustion based on the PM combustion 
technology. Durst [16] proposed two structural versions for PM en¬ 
gine design. In order to prove the feasibility of the PM engine, they 
modified a single-cylinder, air cooled diesel engine to incorporate a 
PM in the cylinder head and operated it as a PM engine .The engine 
was operated for several hours without any damage to the PM 
material or showing any uncontrolled operations. The experiments 
showed very low emission level for the test engine and much lower 
combustion noise due to the reduction in the pressure peak com¬ 
pared to conventional DI engines. However, no data for perfor¬ 
mance of PM engines are available in the open literature. 

Referring to the design of Durst [16], this paper presents a two- 
zone model to investigate the combustion process of the PM en¬ 
gine. A computer program was developed and coupled with the 
chemical kinetics package CHEMKIN III. In this way a skeletal ki¬ 
netic mechanism for iso-octane oxidation including 38 species 
and 69 elementary reactions was used for the chemistry simula¬ 
tion, which could predict satisfactorily ignition timing, burn rate 
and the emissions of HC, CO and NO* for HCCI engine [17]. 



Fig. 1. Principle of the PM-engine cycle. (Left) permanent contact PM engine, and 
(Right) period contact PM engine. 


heat released from the combustion process is partly deposited in 
the PM volume, which could be used to vaporize fuel and preheat 
the mixture in the next cycle. The cylinder and the PM chamber 
keep separated during the intake, compression and exhaust 
strokes, which make this type of PM engines operate as a conven¬ 
tional DI engine in these strokes. 

The main difference between the two type PM engines is 
whether the compressed air contacts with fuel vapors during its 
evaporation, which decides the ignition timing. 

The formation of homogeneous mixture and 3D-thermal self¬ 
ignition in the PM volume create a realizable condition for homo¬ 
geneous combustion, which is almost independent of the engine 
load. Furthermore, the heat recuperation in the PM may be used 
for heating up the compressed air and controlling the combustion 
temperature level. 

2.2. Computational model 


2. Mathematical model 

2.1. Working principle of PM engine 

The PM engine is here defined as an internal combustion engine 
with a highly porous medium (PM) chamber mounted on the cyl¬ 
inder head (Fig. 1). Same components as in conventional internal 
combustion engines are still used in the PM engine, except the fuel 
injection system that is required to be modified and optimized. 

In the “permanent contact” PM engine the PM chamber always 
contacts with the cylinder and its operation is the same with a tra¬ 
ditional diesel engine. In the “period contact” PM engine, the PM 
chamber is equipped with a valve permitting a periodic contact be¬ 
tween the PM chamber and the combustion chamber. At the end of 
expansion stroke the valve is closed and fuel is injected into the PM 
volume, the process continues through exhaust, intake and com¬ 
pression strokes. The large specific surface area of the PM and 
excellent heat transfer in the PM volume as well as a long time 
available for fuel vaporization make the liquid fuel vaporized com¬ 
pletely. At the end of compression (near the TDC), the valve of the 
PM chamber opens and the compressed air flows from the cylinder 
into the hot PM volume containing fuel vapors, generating a highly 
turbulent flow condition. Very fast fuel-air mixing occurs and the 
resulting mixture is self-ignited across the whole PM volume. The 


In the present model the combustion chamber is divided into 
two zones (Fig. 2) with different temperatures and mass composi¬ 
tions according to the structure of the PM engine: the PM chamber 
zone (zone one) and cylinder zone (zone two). During the combus¬ 
tion process the volume of PM zone keeps constant while the vol¬ 
ume of cylinder zone varies with time (crank angle). The 
thermodynamic properties and mass composition are spatially 
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Fig. 2. The definition of two zones in different engine. (Left) in permanent contact 
PM engine, and (Right) in period contact PM engine. 
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uniform in each zone. Moreover, mass exchanged between two 
zones occurs to maintain the in-cylinder pressure uniform. 

For simplicity, the following assumptions are introduced in the 
formulation: 

(a) The porous medium with high porosity is homogeneous, 
continued and non-catalytic [11,13]. 

For the period contact PM engine: (b) Fuel droplets have vapor¬ 
ized completely before the valve is open due to the excellent heat 
properties of PM volume and the long time available for fuel vapor¬ 
ization [15,16], Consequently the temperature of air, fuel vapor 
and the solid in the PM chamber are uniform as the valve is open, 
(c) When the valve is open, the highly compressed air is assumed 
to be mixed with the fuel vapor immediately to form a homoge¬ 
neous combustible mixture because the pressure in cylinder cham¬ 
ber is much higher than that in the PM chamber. 

For permanent contact PM engine: (d) All fuel droplets have the 
same diameter and are uniformly distributed within the PM vol¬ 
ume when fuel is injected in the PM chamber. Through the vapor¬ 
ization process the total number of fuel droplets decreases while 
the diameter keeps constant, which influences the mass fraction 
of fuel vapor [12]. (e) The heat transfer between the PM solid 
and droplets is negligible for little probability of contact between 
each other. 


2.2.2. Thermodynamics and chemical kinetics model 

In the thermodynamic model the mass of each zone is assumed 
to be constant and species are assumed to obey the ideal gas equa¬ 
tion of state. According to the first law of thermodynamics (FLT), 
the energy equations of each zone can be written as: 

a - dTi dVi 

Qi = m i c v ,, w + p cylw (1) 


The in-cylinder pressure at each time step is calculated using the 
ideal gas relation [8] for all the zones: 



( 


2 


5>,R,T, 

!=1 


j (Y cyl + PMb) 


Substituting Eq. (2) into (1) yields: 



For the permanent contact PM engine, the evaporation and 
chemical reaction take place in PM zone at the same time. The 
mass fraction of fuel vapor is controlled by both chemical kinetics 
and evaporation rate of fuel droplets in PM chamber. Thus, when 
i = 1 and j = ifuel a coefficient Sy = 1 is enrolled where ifuel repre¬ 
sents the species of fuel, and Sy = 0 otherwise. rh f is the evapora¬ 
tion rate of fuel droplets, and is given in model of evaporation. 

For the period contact PM engine Sy = 0, because the evapora¬ 
tion is separated from the chemical reactions and has been finished 
before the valve is open. 

The chemical heat release of zone i is obtained as: 

K dY K 

Qchem,i = tljj ^ Ujj(d)jj ViWjj + Sy • 171 f Vi ) (5) 

J=1 j= 1 

where u jt i is specific internal energy of species j in zone i. Sf has the 
same value as in Eq. (4). 

The overall wall heat loss q /oss/ can be determined by a heat 
transfer model for the HCCI engine, especially. The heat transfer 
in the PM q PM j can be obtained similarly by a heat transfer model 
for the PM. 

Eq. (3) provides two equations for the two zones. The variation 
of the mean temperature in each zone dTJdt can then be calculated 
when parameters m h R t , T u c ph q chem ,u Qpm,u Qioss.i for the gas mixture 
and Vcyi , dV C yildt for the engine are given. Combining with the spe¬ 
cies mass fraction (Eq. (4)) for each zone, an equation set of 2 (K + 1) 
order is built. Giving the initial conditions, variations in tempera¬ 
tures and species mass fractions can be obtained by solving the 
equation set with the DVODE package from CFIEMKIN which is a 
variable-coefficient Ordinary Differential Equation solver, with 
fixed-leading coefficient implementation. 

In the above calculation, the mass transferred between zones is 
not considered. Thus, a special model is needed for computing the 
mass exchange at the end of each time step to redistribute the 
mass of each zone and keep the in-cylinder pressure uniform. 

For HCCI engines, it was found that the time step of 0.5 CAD is 
small enough [9]. In the following calculations, the time step is ta¬ 
ken to be 0.4 CAD before ignition and 0.06 CAD during the combus¬ 
tion process, which is adequate for predicting the temperature 
distribution and evolution. 


dl) 

dt 



Qchem , i + Qloss,i + QPM,i 


T.Lm,R,T, dV, 

i7ii(V cy i + V PM b) dt 



In Eqs. (1)—(3), m it T„ V» c v>i , Ri and Qi are mass, temperature, vol¬ 
ume, specific heat capacity at constant pressure, universal gas con¬ 
stant and heat release in each zone i, respectively. V PMB is volume of 
the PM chamber, and s is the porosity of the PM. The cylinder vol¬ 
ume V cy i and its variation rate dV cy ildt can be described by functions 
of crank angle and cylinder size [18]. For unit mass of the mixture in 
the single zone i: q chem , i is heat release of chemical reaction, q /oss , i is 
heat loss through the cylinder wall into the surroundings, q PM ,i is 
the total heat transfer in the PM between gas and liquid as well 
as gas and solid, where dVf \dt = 0 for 2 = 1, dVJdt = dV cy ildt and 
Qpm,\ = 0 for i = 2 . 

The chemical heat release in each zone is obtained by the fol¬ 
lowing method. 

For each species j in zone 2, a set of differential equations for 
species mass fraction Y jti is solved using the CHEMKIN III libraries 
to obtain the change in mixture composition [8,9]: 

dY 

d( ! COliViWjj I (4) 


2.2.2. Fteat transfer model for cylinder wall 

It is well known that modeling of wall heat losses in HCCI com¬ 
bustion processes is a very challenging task. The overall wall heat 
loss is defined as: 

d J v = —hA(T - T w ) (6) 

where h , A and T w are heat rate coefficient, combustion chamber 
area and wall temperature, respectively. 

According to the experimental results reported for the gasoline 
HCCI process, the improved Woschni model of Chang [20] showed 
very good performance for the prediction of heat transfer in HCCI 
engine. The heat transfer rate coefficient is written as: 

h = \29.8B~° 2 P° a r 0 73 co 08 with: 

where B is the cylinder bore, T is mean gas temperature, Ci and C 2 
are constants in the Woschni heat transfer model and co is mean gas 
velocity defined as function of average velocity of piston 2), cylinder 
volume V d , referenced temperature T r , pressure p r , volume V n and 
motoring pressure p mot . 


where djjj is the molar rate of production calculated by CHEMKIN 
[19]. Vi is the average specific volume of zone 2. W jti is molecular 
weight of species j in zone 2. 


2.2.3. Modeling of fuel evaporation 

223A. For permanent contact PM engine. Since modern nozzles are 
able to produce a fuel spray with a mean diameter in the order of 
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10 pm, the initial diameter 100 jam is considered in this study. 
Therefore, the average size of the fuel droplets is small enough 
compared to the average pore size of the PM. As assumed, the con¬ 
tact heat transfer between liquid and solid is neglected. Further¬ 
more, the radiative heat transfer of the solid to the droplets is 
not taken into account as the temperature of gas is higher than that 
of the solid during the fuel injection. The heat transfer between gas 
and solid is dominant, the heat transfer between gas and liquid are 
also considered. 

When the temperature of liquid is lower than its saturation 
temperature, the energy equation for liquid phase [12] is written 
as: 


epfCi 


dh 

dt 


y g iTg - t,) 



where y gl is the volumetric heat transfer coefficient between gas 
and liquid. T t is the temperature of liquid droplet. p f = p,n(7rd 3 /6) 
is the density of fuel, d and n is the mean diameter and the number 
density of liquid droplets, respectively, pi is the material density of 
liquid. 

The volumetric heat transfer coefficient between gas and liquid 
per unit volume is defined as y gl = n ■ s • h gt • S d [12], where h g[ is 
convective heat transfer coefficient between liquid and gas and 
S d is the specific surface of a liquid droplet. The number density 
of liquid droplet n is given as: n = , ™ Uq r uid w —, where m/ I0Uld is 

the mass of liquid fuel contained in the space of porous medium. 
The initial value of n is determined by the total mass of fuel in¬ 
jected per cycle and a given droplet diameter. With the progress 
of evaporation, n decreases until to zero. Because inside the PM 
turbulence has only a weak effect, assuming the Nusselt number 
[12] Nu g i = h g id/l g = 2, and h gi can be calculated. Thus, the volu¬ 
metric heat transfer coefficient between gas and liquid can be de¬ 
scribed as: 


y g i = ^ nndNUgi • X g 

The evaporation rate of fuel droplets m f is given as: 

fh/ = y g i(T g — T sat )/L 


( 8 ) 

( 9 ) 


where y gi is the volumetric heat transfer coefficient between gas 
and liquid, T g is the temperature of gas mixture, T sat is the saturation 
temperature of solvent-fuel mixture and L is the latent heat of 
evaporation. 


2.23.2. For period contact PM engine. For the period contact PM en¬ 
gine, the evaporation process of liquid fuel occurs in PM chamber 
during the valve keeps closed. Thus, the evaporation process is 
decoupled from the cylinder, and it is not necessary to study the 
detailed evaporation process. The most important problem is then 
the thermodynamic properties of gas mixture and solid matrix at 
the time of the valve opening. Here, it is assumed that the temper¬ 
ature of porous solid and fuel vapor at the valve opening are uni¬ 
form because a long time available for evaporation and heat 
transfer within the PM. The temperature inside a fuel droplet is 
also assumed uniform, but varies with time. The stress of the evap¬ 
oration process is put on calculation of the heat transfer between 
gas-solid and gas-liquid in PM chamber since the valve was 
closed, whereas the evaporation rate of liquid droplets is not men¬ 
tioned here. The heat balance equation for the closed PM chamber 
can be formulated as follows: 

p s c s (T sc T open ) + pec(T gc - T open ) 

— ITlfCi (T sat T/) + ITlfL + C evap mf(J open — Tsat) (10) 

where p s and p are the densities of the PM and remained gas in the 
PM chamber, respectively. T sc and T gc are the temperatures of the 


solid and gas in the PM chamber when the valve closed at former 
cycle. T open is the temperature of vapor and solid in PM chamber 
as the valve was open, which can be calculated by Eq. (10). m f is 
the total mass of fuel injected per cycle, which is determined by 
the mass of intake air and a specific equivalence ratio. c s , c, q and 
Cgvap is the specific heat of solid, remained gas, fuel droplet and fuel 
vapor in the PM chamber, respectively. 

2.2.4. Heat transfer modeling in porous medium 
2.2.43. For permanent contact PM engine. For the permanent con¬ 
tact PM engine, the heat transfer associated with liquid droplet 
during evaporation should be considered in the model of heat 
transfer in porous medium. 

The total heat transfer in porous medium q PM ,i can be described 
with the following equation: 

= [Tgs ■ (Ts - Tg) - (1 - S)y gl (T g - T,)-S ■ m r 
L - 5 ■ C evap ■ rh f ■ (T g - T sat )]/p (11) 

where T g and T s are the temperatures of gas mixture in the porous 
space and the solid porous media, respectively. y gs is the volumetric 
heat transfer coefficient between solid and gas. C evap is specific heat 
of fuel vapor and p is the density of gas mixture in porous space. 3 is 
a coefficient of evaporation, 3 = 0 for T t < T sat and 3 = 1 for T t = T sat 
[12]. 

Energy equation for solid phase [21] is: 

Ps c s^ = y gs (T g -T s ) (12) 

In order to get y gs , a volumetric Nusselt number Nu v is defined as 
[12]: 

Nu„ = y gs D 2 /l g (13) 

where D is the average diameter of pore, I g is the thermal conduc¬ 
tivity of gas mixture. The correlation for the volumetric Nusselt 
number Nu v [21] is given by the following correlation: Nu v = CRe m . 
Values of C= 0.187 and m = 1.1 are used, which depend on the pore 
diameter. For diesel engine, the Reynolds number ( Re ) is calculated 
by Re = Bo/v g , where i> is the mean piston speed as the characteris¬ 
tic velocity. B is cylinder bore as the characteristic length, v g is the 
mean kinematic viscosity of gas mixture. 

2.2.42. For period contact PM engine. For the period contact PM en¬ 
gine, the volume of porous medium keeps in contact with the cyl¬ 
inder as the valve to the PM chamber keeps open. In this study a 
model for heat transfer in the PM has been developed, which in¬ 
cludes an additional energy equations for the solid matrix and cal¬ 
culation of heat transfer coefficient between gas and solid matrix. 
Because the liquid fuel has been vaporized completely at the valve 
opening, the heat transfer relate to the liquid phase is neglected. 
The energy equation for solid phase is the same as in the case of 
the period contact PM engine. 

The heat transferred q PMJ to unit mass of the mixture in the PM 
chamber zone can be expressed as: 

Qpm, i = y gs ■ (Ts - Tg)/p (14) 


2.2.5. Model for mass exchange between two zones 

For the both type of PM engines a model for mass transfer between 
two zones is needed when the PM chamber is coupled with the cylin¬ 
der. In the present paper, the mass exchange model is adopted and 
simplified from the model of Komninos and Hountalas [8] based on 
the assumption of uniform pressure throughout both zones. 

At each time step, since the pressure is uniform throughout the 
cylinder, the pressure at each time step is calculated using the ideal 
gas relation for zone i: 
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Pcyi = mjRiTi/Vi, i = l,2 (15) 

Solving Eq. (15) for the mass of each zone and summing for all the 
zones: 


m cyl 


x> 

i=i 


P<yi^2i =i 


A 

RJi 



Solving Eq. (16) for the uniform cylinder pressure: 


2 


Pcyl = J2 m > = 
1=1 


m cy i 

Ai=l Rj i 



where the temperature of each zone T, is calculated by Eq. (3), 14 m, 
and Ri are volume, mass and universal gas constant of each zone. 

At each time step the mass of each zone is calculated using the 
equation of state for ideal gas: 


mi = 


Pcyl^i 


i = 1,2 



The mass change of each zone during the time step is therefore: 

A m i = m i ~mP reviousCA , i = 1,2 (19) 


Calculations in each time step are divided into two parts. First, the 
change of composition of each zone is computed using thermody¬ 
namics and chemical kinetics with the assumption of no mass 
transfer between two zones. Next, at the end of the step, the sub¬ 
routine of mass exchange is called to compute the mass changes 
to keep the pressure uniform throughout the cylinder, and then re¬ 
new associated parameters. Therefore, the change of composition of 
each zone is the result of combustion and mass transfer from an¬ 
other zone. 

The mass flowing from one zone to another has the temperature 
and chemical species composition of the zone from which it origi¬ 
nates. Therefore, the exchange of mass will result the transfer of 
species and enthalpy between different zones. 


2.3. Modeling of chemical kinetics 


Combustion is described using a skeletal chemical kinetics 
mechanism for iso-octane oxidation consisting of 38 species and 
69 reactions [17]. For each species, the rate of production is calcu¬ 
lated and the differential equations of species mass fraction are 
solved using the chemical kinetics package CHEMKIN III to obtain 
the changes in mixture composition in each zone. Comparisons 
with various experimental data, including shock tube, rapid com¬ 
pression machine, jet stirred reactor and HCCI engine, have indi¬ 
cated very good performance of this mechanism over wide 
ranges of temperature, pressure and equivalence ratio, especially 
at high pressure and lean equivalence ratio conditions. Application 
of the skeletal mechanism to a single-zone model of HCCI engine, 
yielded results substantially identical with those from the detailed 
mechanism developed by Jia and Xie [17], and furthermore, the 
computing time was reduced greatly. 

2.4. Engine specifications and operating conditions 


Following the study by Durst [16], the PM engine simulated in 
this paper is based on a Cummins B-series engine, which is a typ¬ 
ical medium-duty diesel engine. Dec and co-worker [7,22] have 
conducted systematical experimental investigation on the engine, 
which was converted for operation in the HCCI mode. The detailed 
Specifications of the engine is showed in Table 1. Following the 
suggestion of Sjoberg and Dec [7,22], an effective compression ra¬ 
tio of 16.8 is used for the calculation. For PM engine, the existence 
of the PM chamber enlarges the total space volume of combustion 
chamber. Sequentially, the mass of intake air for the PM engine is 


Table 1 

Engine specifications. 


Displacement(single-cylinder) 

0.981 1 

Bore 

102 mm 

Stroke 

120 mm 

Connecting rod length 

192 mm 

Nominal geometric compression ratio 

18:1 

Fuel 

Iso-octane 

Cylinder wall temperature 

384.09 K 


0 CA is taken to be TDC intake, so TDC combustion is 360°. 


larger and the actual compression ratio is smaller than those of the 
original engine if other parameters keep unchanged. 

The baseline operation condition is specified according to the 
original engine: the intake pressure is maintained at 100 kPa 
throughout the study, simulating naturally aspirated operation; 
the intake temperature of air is 320 K, much lower than that in 
HCCI (400 I<) [7,8]; an engine speed of 1200 rpm is used for all 
operations; the excess air ratio is 1.6; the injection occurs at 352 
crank angle for the permanent contact PM engine. Parameters rel¬ 
evant to the PM chamber for both types PM engine are: the initial 
temperature of the PM volume is 1100 K, which is approximately 
the average in-cylinder temperature of a common diesel engine; 
The cylindrical PM chamber mounted on the cylinder head is 
40 mm in diameter and 50 mm high; the material of porous cera¬ 
mic foams is silicon carbide with a porosity of s = 0.87, density of 
Ppm = 510 kg/m 3 , porous diameter of D = 1.52 mm and specific heat 
of C PM = 824.7 J/kgK [21]. For the convenience of comparison, all 
the operating parameters can change separately in the calculation. 

3. Results and discussion 

It is known that the auto-ignition temperature of hydrocarbon 
fuels used in engines has often been determined according to the 
sudden increase in the in-cylinder (including both zones) gas pres¬ 
sure or temperature [7]. As key parameters for controlling the 
auto-ignition process the change in the in-cylinder gas pressure 
or temperature are investigated with emphasis in this paper. 

3.1. Influence of the initial temperature 

Fig. 3 shows effect of the intake temperature on the in-cylinder 
temperature and pressure in permanent contact PM engine. It can 
be seen that the crank angle of the auto-ignition in the cylinder 
zone (zone two) varies significantly with the initial temperature. 
However, the process of combustion in the PM zone (zone one) 
shows almost no change with the change in the initial tempera¬ 
ture. That is because the hot PM preheats the gas in PM zone, the 
temperature of gas mixture in the PM zone before combustion is 
mainly affected by the temperature of PM, which is relatively stea¬ 
dy for each cycle. 

If the initial temperature is below a certain threshold, the auto¬ 
ignition can occur only so late or even not occur that a sharp rise in 
the pressure curves cannot be observed. The simulation has shown 
that the minimum initial temperature required for an auto-ignited 
combustion in the cylinder zone with the present engine configu¬ 
ration and operating conditions should be above 300 K. With 
increasing the initial temperature, an obvious two-stage ignition 
process including both low and high temperature kinetics is dis¬ 
played in Fig. 3. Due to the combustion in the PM zone, the temper¬ 
ature of mixture in the PM zone increases and mass flows from the 
PM zone into the cylinder zone as computed in the mass transfer 
model. Subsequently, the mixture in the cylinder zone is ignited, 
low temperature reactions occur and follow by high temperature 
reactions after a long time delay. With the change in temperature, 
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Fig. 3. Effect of the initial temperature on cylinder temperature (a) and pressure (b). 


pressure also increases in two stages. The ignition process of the 
PM engine in the cylinder zone is similar to the HCCI engine. 

For a HCCI engine, a fully premixed fueling system is necessary 
to supply a homogeneous combustible mixture and the initial tem¬ 
perature must be high enough to enable an auto-ignition at the end 
of compression stroke. In contrary, for the PM engine fuel can be 
just injected directly into the PM chamber with a lower initial tem¬ 
perature (e.g. 320 K) and the evaporation process occurs in the PM 
volume and is separated from the cylinder. This is an outstanding 
advantage of the PM engine. The initial temperature of 320 K was 
used in the following simulations, unless otherwise stated. 

3.2. Influence of the compression ratio 

Fig. 4 shows effect of the compression ratio on the in-cylinder 
temperature in permanent contact PM engine. It can be seen that 
auto-ignition in the cylinder zone could occur if compression ratio 
is higher than 13. With increasing compression ratio the ignition 
point moves ahead and the maximum temperature in the cylinder 
increases. The calculation indicates that the influence of compres¬ 
sion ratio on PM engine is not as strong as on the common HCCI 



Fig. 4. Influence of the compression ratio on cylinder temperature. 


engine. As the combustion starts in the PM zone, mass and heat 
are transferred from the PM zone to the cylinder zone, which 
makes ignition and combustion of mixture in the cylinder more 
easily. In order to attain a reasonable compromise between the 
power output and emissions of NO* and HC, it is still important 
to choose a proper compression ratio to keep the ignition timing 
nearly at the TDC. 

3.3. Influence of the initial pressure 

Fig. 5 shows the effect of the initial pressure on the in-cylinder 
temperature in permanent contact PM engine. For high initial pres¬ 
sure, auto-ignition in the both zones move ahead, especially in the 
cylinder zone, i.e. the ignition delay between the two zones gets 
shorter. That is because increasing initial pressure enlarges corre¬ 
spondingly the fuel mass, which results in more combustion heat 
released in the PM zone and much earlier ignition in the cylinder 
zone. 

In Fig. 5, it is shown that with increasing initial pressure the 
maximum cylinder pressure increases evidently while the maxi¬ 
mum temperature almost shows no change. This means that the 
thermal efficiency and NO* emission, being effected greatly by 
peak temperature, will remain at a relatively steady level when 
the initial pressure increases. Moreover, the voids in the PM results 
a lower pressure in the PM engine than that in the HCCI engine, 
which is beneficial for preventing engine knock. It can therefore 
be anticipated that boosting intake pressure is an attractive meth¬ 
od in the PM engine for higher power output. 

3.4. Influence of the excess air ratio 

Fig. 6 shows effect of the excess air ratio 2 on the in-cylinder 
temperature and pressure in permanent contact PM engine. For 
2 = 1.0, the peak temperature and pressure in the two zones are 
fairly high and the ignition delay between two zones reaches the 
minimum. With increasing the excess air ratio, the curves of tem¬ 
perature and pressure trend smooth and the auto-ignition in the 



Fig. 5. Influence of the initial pressure on cylinder temperature (a) and pressure (b). 
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Fig. 6. Effect of the excess air ratio on cylinder temperature (a) and pressure (b). 


cylinder occurs very late, until it cannot ignition for 2 = 2.0. Calcu¬ 
lation indicates that the combustion efficiency increases with 
increasing the excess air ratio, because a leaner fuel-air mixture 
could combust more completely. Furthermore, it is easy to under¬ 
stand that the exhaust concentrations of NO, C0 2 and CO, as func¬ 
tions of the excess air ratio, will also decrease. 

In conclusion, a proper excess air ratio must be used to ensure 
both the occurring of ignition and low emissions for the PM engine. 

3.5. Influence of the engine speed 

Fig. 7 shows effect of the engine speed on the in-cylinder tem¬ 
perature in permanent contact PM engine. It is showed that the 
auto-ignition in the PM zone delays about 5 CA and the interval be¬ 
tween ignition timings in the two zones gets longer with engine 
speed increasing from 800 to 2000 rpm. An auto-ignition in the 
cylinder cannot happen when engine speed reaches 2000 rpm. 
With increasing engine speed, less time is available for the auto¬ 
ignition reactions to take place near the TDC, thus the combustion 
of mixtures becomes more incompletely. As a result of this, the 
peak combustion temperature in the PM zone decreases with in¬ 
creased engine speed for the same excess air ratio. It can be noted 
that for a high speed PM engine one should take some measures in 
order to realize a complete combustion near TDC, such as boosting 
pressure or increasing compression ratio. 

3.6. Influence of the initial temperature of porous medium 

Fig. 8 shows effect of the initial temperature of the porous med¬ 
ium insert on the in-cylinder temperature in permanent contact 
PM engine. It is clear from Fig. 8 that the PM temperature is a dom¬ 
inating parameter determining the auto-ignition in the PM zone, 
combustion cannot occur in the PM zone or in the cylinder zone 
when the initial PM temperature is lower than 900 K. At 1000 K, 
the combustion occurs in the PM zone, however, ignition in the cyl¬ 
inder zone still could not be realized. Only with farther increase in 



Fig. 8. Effect of the PM temperature on cylinder temperature. 

temperature auto-ignition in both zones can succeed. The calcula¬ 
tion indicates that the higher the initial PM temperature is, the eas¬ 
ier the auto-ignition takes place, at the same time more NO x 
emission will be produced. In conclusion, the initial PM tempera¬ 
ture is an important factor determining the possibility of ignition, 
the crank angle of auto-ignition for each zones and the average in¬ 
cylinder temperature, which is an important object for PM engine 
research. 

3.7. Mass exchange between zones 

Fig. 9 shows variations of the masses contained inside the cylin¬ 
der zone (zone 2) and PM zone (zone 1) in permanent contact PM 
engine. The total mass of mixture gas increases during fuel evapo¬ 
ration process, which is shown in detail. In Fig. 9, it is shown that 
the mass of mixture is mainly distributed in cylinder zone; during 
compress stroke compressed air flows from the cylinder zone into 
the PM zone and opposite variety happens during expansion 
stroke. As combustion occurring in the PM zone near TDC, the tem¬ 
perature and pressure in the PM zone rise sharply and mass flows 
into the cylinder zone, when the combustion in cylinder happens, 
mass flows from the cylinder zone to the PM zone. That means that 




Fig. 7. Effect of the engine speed on cylinder temperature. Fig. 9. Mass distributions of two zones. 
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Fig. 10. Effect of the valve opening time 
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the auto-ignition in both zones is the crucial factor for the mass ex¬ 
change between two zones. Furthermore, the combustion is con¬ 
trolled by the chemical kinetics model, which means the more 
fuel mass result the more heat it released, and the calculation also 
shows that the heat released has same trend with the mass 
distributing. 

3.8. Influence of the valve opening time in the period contact PM 
engine 

For the period contact PM engine, The Influence of the varied 
parameters (e.g. the compression ratio, engine speed and excess 
air ratio, etc.) on the temperature of the in-cylinder temperature 
are approximate to that in permanent contact PM engine, which 
are not discussed in detail here. 

Fig. 10a shows effect of the valve open timing on the in-cylinder 
temperature in the period contact PM engine. The excess air ratio is 
2.0 and the other parameters are same with those in permanent con¬ 
tact PM engine. The maximum temperature is much lower than that 
in permanent contact PM engine, for a larger excess air ratio is used. 
It is shown that the combustion in the PM zone follow the valve open 
and the time delay is about 4 CA. With the increase of valve opening 
time, the ignition delay between two zones gets shorter. 

It is worth to note that, in Fig. 10b there is a sharply drop in pres¬ 
sure when the valve opens. Because the pressure in the PM cham¬ 
ber, which remains its value at the end of former cycle, is much 
lower than that in the cylinder before the valve opening. The low 
pressure continues until the combustion in the PM zone happens, 
and then the pressure increases. The wave of pressure in cylinder 
might cause shake of the engine, which should be avoided. More¬ 
over, the combustion in cylinder zone has not occurred as the valve 
open at 340 CA, because the air in cylinder has not got enough com¬ 
pression and its temperature is not high enough to cause auto¬ 
ignition. 

4. Conclusions 

A two-zone model for the new concept internal combustion en¬ 
gine based on the technology of combustion in porous medium is 
presented in this paper. The main attention is focused on its poten¬ 
tial for the realization of a homogeneous combustion process un¬ 
der variable engine operational conditions. The results obtained 
can be summarized as follows: 

1. The PM engine is fueled by direct injection of liquid fuel into the 
PM chamber with a low initial temperature (300 I<) and the 
evaporation process occurring in the PM volume is separated 
from the cylinder, for which a premixed fueling system and suf¬ 
ficiently high initial temperature are not need as in the case of 
HCCI engine. 


2. The compression ratio of 13 is great enough for the auto-igni- 
tion in the PM zone, which is not so crucial like in the HCCI 
engine, because the porous medium with high temperature pre¬ 
heats the combustible mixture. Besides, the combustion heat 
released in the PM enhances the ignition in the cylinder. This 
means that the PM engine could have a wider operation range 
than a HCCI engine. 

3. The auto-ignition timing in the PM zone is dependent mainly on 
the initial PM temperature, ignition can easily occur as the PM 
temperature exceeds 1000 K depending on operating condi¬ 
tions. Moreover, the combustion process in the PM zone is rel¬ 
atively steady, which can be controlled by adjusting the initial 
PM temperature. 

4. For the same intake pressure of 1 atm, the pressure in the PM 
engine is relatively lower than that of a traditional engine due 
to the voids in PM, boosting in the PM engine may shorten 
the ignition delay and enlarge the peak pressure with little 
effects on the peak temperature, which is an attractive method 
to improve engine output. 

5. The average temperature in the PM zone is higher than that in 
the cylinder, resulting in slightly more emissions of NO for 
existing excess air ratio of 1.6, which is still a problem for PM 
engine. In the future research more attention should be paid 
to lean combustion in the PM engine. 
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